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A phylogenetically distinct Middle East respiratory 
syndrome coronavirus detected in a dromedary calf 
from a closed dairy herd in Dubai with rising 
seroprevalence with age 
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Alan KL Tsang’, Nissy Annie Georgy Patteril’, Honglin Chen*®”*, Shyna K Elizabeth’, Kwok-Yung Yuen””””, 
Sunitha Joseph’, Ningshao Xia’, Renate Wernery’, Susanna KP Lau”®””” and Patrick CY Woo”®””* 


Middle East respiratory syndrome coronavirus (MERS-CoV) was detected by monoclonal antibody-based nucleocapsid protein-capture 
enzyme-linked immunosorbent assay (ELISA), RNA detection, and viral culture from the nasal sample of a 1-month-old dromedary calf 
in Dubai with sudden death. Whole genome phylogeny showed that this MERS-CoV strain did not cluster with the other MERS-CoV 
strains from Dubai that we reported recently. Instead, it formed a unique branch more closely related to other MERS-CoV strains from 
patients in Qatar and Hafr-Al-Batin in Saudi Arabia, as well as the MERS-CoV strains from patients in the recent Korean outbreak, in 
which the index patient acquired the infection during travel in the eastern part of the Arabian Peninsula. Non-synonymous mutations, 
resulting in 11 unique amino acid differences, were observed between the MERS-CoV genome from the present study and all the other 
available MERS-CoV genomes. Among these 11 unique amino acid differences, four were found in ORF Lab, three were found in the S1 
domain of the spike protein, and one each was found in the proteins encoded by ORF4b, ORF5, envelope gene, and ORF8. MERS-CoV 
detection for all other 254 dromedaries in this closed dairy herd was negative by nucleocapsid protein-capture ELISA and RNA 
detection. MERS-CoV IgG sero-positivity gradually increased in dromedary calves with increasing age, with positivity rates of 75% at 
zero to three months, 79% at four months, 89% at five to six months, and 90% at seven to twelve months. The development of a rapid 
antigen detection kit for instantaneous diagnosis is warranted. 
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INTRODUCTION 

Since its first appearance in 2012, the Middle East respiratory syn- 
drome (MERS) has affected more than 25 countries in four continents, 
with more than 1300 cases and a frightening fatality rate of greater than 
30%.’ A novel lineage C betacoronavirus the MERS coronavirus 
(MERS-CoV), has been confirmed to be the etiological agent.” The 
subsequent detection of MERS-CoV and its antibodies in dromedaries 
in various countries in the Middle East has led to the hypothesis that 
these animals may be the reservoir for MERS-CoV.* '* Furthermore, 
the discovery of other closely related lineage C betacoronaviruses both 
before and after the beginning of the MERS epidemic in various bat 
species and hedgehogs, '*~'® as well as the demonstration of the ability 
of the spike protein of Tylonecteris bat CoV HKU4 to bind to dipepti- 
dyl peptidase 4 (DPP4), the receptor of MERS-CoV, have also sug- 
gested that these animals may be the hosts for the ancestor of MERS- 
Cov, 19-22 


Although MERS-CoV can be cultured and antibody detection 
methods are available, laboratory diagnosis and epidemiology studies 
of MERS-CoV infections are mainly achieved by real-time quantitat- 
ive reverse transcription polymerase chain reaction (real-time qRT- 
PCR), using the RNA-dependent RNA polymerase, the region 
upstream of the envelope gene or the nucleocapsid (N) gene of 
MERS-CoV as the target.” However, the cost of real-time qRT- 
PCR is still high and such expertise may not be available, particularly 
in some clinical microbiology laboratories in the Middle East where 
most of the cases were identified. Because enzyme-linked immuno- 
sorbent assay (ELISA) for antigen detection offers an inexpensive and 
user-friendly approach for laboratory diagnosis and epidemiological 
studies of respiratory viral infections, including severe acute respir- 
atory syndrome,~*”” we recently developed a monoclonal antibody- 
based capture ELISA for the detection of MERS-CoV N protein.*° In 
this study, we describe the application of N protein ELISA, real-time 
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qRT-PCR, isothermal amplification, virus isolation, and antibody 
detection in a cross-sectional epidemiological study of a dromedary 
dairy herd in Dubai after sudden death and the detection of MERS- 
CoV from a dromedary calf within the herd. 


MATERIALS AND METHODS 


Animals and samples collection 

The 1-month-old female dromedary calf with sudden death belonged 
to a closed dairy herd reared for milk production, which roamed a 
small desert area in Dubai, the United Arab Emirates. In addition to 
the calf that succumbed, samples were collected from the remaining 
254 dromedaries in the herd, comprising 133 adult female dromedar- 
ies and 121 dromedary calves. There were no adult male dromedaries 
in the herd and mating was human-assisted. 


Capture ELISA for detection of MERS-CoV N protein 

MERS-CoV N protein ELISA on nasal swabs was performed as prev- 
iously described.*° Briefly, microplates (Sigma-Aldrich, St. Louis, 
MO, USA) were pre-coated with the first anti- MERS-CoV-recombin- 
ant N protein monoclonal antibody 1F6 and incubated at 37 ‘C over- 
night with a blocking reagent (phosphate-buffered saline with 2% 
sucrose, 0.2% casein-Na, and 2% gelatin). Fifty microliters of viral 
lysis buffer was added to each coated well and 50 uL of each nasal 
sample was then added to the wells in duplicate. The plate was shaken 
for 2 min and incubated at 37 ‘C for 30 min. After the wells were 
washed five times, 100 uL of the second anti- MERS-CoV-recombinant 
N protein monoclonal antibody 7C4 conjugated with horseradish 
peroxidase was added and the plate was incubated at 37 ‘C for 30 
min. After five washes, detection was carried out by adding 100 uL 
of tetramethylbenzidine per well, incubating the samples for 10 min, 
and adding 50 uL of 0.2 M H,SO,. The optical density (OD) 450/630 
nm was measured with a microplate reader. 


Real-time qRT-PCR screening assay 

MERS-CoV nucleic acid detection on the nasal swabs and milk was 
performed using a real-time qRT-PCR assay targeting the region 
upstream of the envelope gene and by isothermal amplification using 
a Genie instrument (OptiGene Limited, Horsham, UK). For real-time 
qRT-PCR assays, a 25-uL reaction was set up containing 5 ul of RNA, 
12.5 uL of 2X reaction buffer provided with the Superscript II one- 
step RT-PCR system with Invitrogen Platinum Taq Polymerase 
(Thermo Fisher Scientific, Inc., Waltham, MA, USA); containing 0.4 
mM of each dNTP and 3.2 mM magnesium sulfate), 1 uL of reverse 
transcriptase/Taq mixture from the kit, 0.4 uL ofa 50 mM magnesium 
sulfate solution (Thermo Fisher Scientific, Inc.), 1 ug of non-acety- 
lated bovine serum albumin (Sigma-Aldrich), 400 nM concentrations 
of primer upE-Fwd (GCA ACG CGC GAT TCA GTT) and primer 
upE-Rev (GCC TCT ACA CGG GAC CCA TA), as well as a 200 nM 
concentration of probe upE-Prb (6-carboxyfluorescein (FAM))-CTC 
TTC ACA TAA TCG CCC CGA GCT CG-6-carboxy-N,N,N,N’ 
-tetramethylrhodamine (TAMRA)). The thermal cycling protocol 
was 55 C for 20 min, followed by 95 C for 3 min and then 45 cycles 
of 95 C for 15 s and 58 C for 30 s. 


Viral culture 

MERS-CoV isolates from nasal swabs were grown on African green 
monkey kidney (Vero) cells.” Samples were inoculated in Vero 
cells at a multiplicity of infection of 1.0 for 1 h. Unattached viruses 
were removed by washing the cells twice in serum-free Gibco 
minimum essential medium (Thermo Fisher Scientific, Inc.). The 
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monolayer cells were maintained in minimum essential medium with 
1% fetal calf serum. All infected cell lines were incubated at 37 C for 
5 days. Cytopathic effects were examined on days 1, 3, and 5 by 
inverted light microscopy. 


MERS-CoV IgG antibody detection 

All serum samples were tested for MERS-CoV IgG antibodies using 
Euroimmun ELISA (Euroimmun AG, Lübeck, Germany). Samples 
were applied at a 1:101 dilution and incubated in the wells of micro- 
titer strips coated with the purified spike protein subdomain S1 anti- 
gen of MERS-CoV. A ratio of >1.0 was considered to be positive. 


Complete genome sequencing and analysis 

The complete genome of the MERS-CoV isolate was amplified and 
sequenced as described previously.°* °° RNA extracted from a nasal 
specimen was used as the template. Sequences were assembled and 
edited to produce final sequences. For phylogenetic analysis, a max- 
imum likelihood phylogenetic tree with 1000 bootstraps was con- 
structed using PhyML v3.0, based on complete concatenated coding 
sequences of 68 available human MERS-CoV and 18 camel MERS- 
CoV genomes. A general time-reversible model of nucleotide substi- 
tution with estimated base frequencies, a proportion of invariant sites, 
and a gamma distribution of rates across sites was applied to conduct 
maximum likelihood analysis. 


Unique mutations in the genome 

The nucleotide sequences of the genomes and the deduced amino acid 
sequences of the MERS-CoV sequenced in this study were compared 
to 68 available human MERS-CoV and 18 camel MERS-CoV gen- 
omes. The nucleotide and amino acid sequences of each coding 
sequence were aligned using MAFFT version 7.°” Unique mutations 
resulting in amino acid substitutions at the coding sequences were 
revealed. 


Nucleotide sequence accession numbers 

The nucleotide sequence of the genome of the camel MERS-CoV in 
this study has been submitted to the GenBank sequence database 
under accession NO KT751244. 
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Figure 1 Sero-positivity of the 254 dromedaries in the herd tested in the 
present study. 


RESULTS 


Necropsy, capture ELISA for detection of the MERS-CoV N protein, 
real-time qRT-PCR screening assay and viral culture 

Necropsy of the dromedary calf with sudden death showed that the 
cause of death was cryptosporidiosis and colisepticemia. The nasal 
sample, but not trachea, lung, and tonsil samples collected from the 
dromedary calf were positive for the N protein (OD 3.712, OD of 
negative control 0.065), real-time qRT-PCR, and virus isolation. In 
contrast, the nasal swabs of the other 254 dromedaries as well as the 
milk produced by the 133 female adult dromedaries were all negative 
for the N protein, real-time qRT-PCR, and virus isolation. 


MERS-CoV IgG antibody detection 
MERS-CoV IgG was positive in the sera of 132 (99%) adult drom- 
edaries and 102 (84%) dromedary calves. This IgG positivity showed a 
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gradual increase in dromedary calves with increasing age, with posi- 
tivity rates of 75% in dromedary calves of age 0-3 months, 79% at 
4 months, 89% at 5-6 months, and 90% at 7-12 months (Figure 1). 


Complete genome sequencing and phylogenetic analysis 

Complete genome sequencing of the MERS-CoV isolated from the 
nasal swab of the dromedary calf with sudden death, as well as phylo- 
genetic analysis, showed that this MERS-CoV strain did not cluster 
with the other MERS-CoV strains from Dubai that we reported 
recently (Figure 2).°* Instead, it formed a unique branch more closely 
related to other MERS-CoV strains from patients in Qatar and Hafr- 
Al-Batin in Saudi Arabia. Notably, the MERS-CoV strains from 
patients in the recent Korean outbreak, in which the index patient 
acquired infection during travel in the eastern part of the Arabian 
Peninsula, also belonged to this cluster (Figure 2). 






Saudi Arabia 





Figure 2 (A) Maximum likelihood phylogeny based on complete concatenated coding sequences of 68 available human MERS-CoV and 18 camel MERS-CoV 
genomes. A general time-reversible model of nucleotide substitution with estimated base frequencies, a proportion of invariant sites, and a gamma distribution of rates 
across sites was applied to conduct maximum likelihood analysis. Bootstrap values (1000 replicates) are shown next to the branches. A scale bar indicates the number 
of nucleotide substitutions per site. The MERS-CoV strain from the present study is shown in blue, those from eastern and western Saudi Arabia are in magenta and 
green, respectively, other strains from the United Arab Emirates are in orange, and those from the recent Korean outbreak are in purple. (B) Map of the Arabian 
Peninsula, with the colors of different geographical regions corresponding to those in panel A. 
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Figure 3 Map of unique mutations in the camel MERS-CoV genome, with the nucleotide positions shown on the scale above it. Blue arrows indicate unique mutations 
observed in the MERS-CoV genome. Amino acid changes are shown above the arrows. 


Unique mutations in the genome 

Using the available complete genomes of MERS-CoV as references, 
non-synonymous mutations, resulting in 11 unique amino acid dif- 
ferences, were observed between the MERS-CoV genome from the 
present study and all the other available MERS-CoV genomes 
(Figure 3). Among these 11 unique amino acid differences, four were 
found in ORF lab (E1059D, P1951L, P2269S, and M61171), three were 
found in the S1 domain of the spike protein (S28P, S459T, and P710S), 
and one each was found the proteins encoded by ORF4b, ORF5, the 
envelope gene and ORFS. 


DISCUSSION 

We report the first cross-sectional epidemiology study of MERS-CoV 
in a closed dromedary dairy herd using antigen, nucleic acid, and 
antibody detection and virus isolation. Similarly to results from pre- 
vious reports, MERS-CoV was detected in a dromedary calf.*”*° This 
dromedary calf with MERS-CoV did not show any respiratory clinical 
signs and had died of cryptosporidiosis and colisepticemia. This is in 
line with our observation that other dromedary calves with sudden 
death and MERS-CoV detected have also died of various diseases, such 
as salmonellosis, clostridial enterotoxemia, and isoporiosis instead of 
pneumonia, and MERS-CoV was detected in the animals only incid- 
entally (unpublished data). Notably, the MERS-CoV strain isolated in 
the present study clustered with other patient strains from east Saudi 
Arabia and Qatar, including the one from a Korean patient who 
acquired infection in this area. This is in contrast to the MERS-CoV 
strains reported in our previous study in Dubai,*® which formed 
another cluster not closely related to the present strain. This indicated 
that more than one strain of MERS-CoV is circulating in Dubai. 
Notably, among the three amino acid changes located at the S1 
domain of the spike protein as a result of non-synonymous mutations, 
none are located at the receptor-binding subdomain, suggesting that 
the camel MERS-CoV is still likely to bind human receptor DPP4. The 
significance of these mutations is unclear and needs to be investigated 
further. The dromedaries in this study belonged to a closed herd that 
did not have contact with other herds because they were browsing in a 
confined area. It is believed that the virus could be circulating in the 
herd, which was probably infected in Saudi Arabia several years ago 
when the dams were calves. Alternatively, the deceased calf may have 
been infected from an unknown external source. 

The present data showed increasing sero-positivity with advancing 
age in dromedary calves. It has been reported that the prevalence of 
MERS-CoV antibodies is significantly higher in dromedary adults 
than in calves.*”*"'? In this study, we further demonstrated that the 
seroprevalence in dromedary calves increased gradually as their ages 
increased, with a sero-positivity of almost 100% in adult dromedaries. 
This is in line with the observation that MERS-CoV infections occur 
almost exclusively in dromedary calves after their maternal antibodies 
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have disappeared. In fact, we have tested samples from the nose, tra- 
chea, tonsil, lung, uterus, and placenta of more than 100 adult female 
dromedaries dispatched to our veterinary laboratory for necropsy, and 
all have shown negative results (unpublished data). 

The high specificity of the N protein ELISA is confirmed in this field 
study. Recently, we reported the development of a monoclonal anti- 
body-based capture ELISA for MERS-CoV N protein and evaluated its 
specificity using nasopharyngeal specimens from patients with various 
common respiratory virus infections, including the four human CoVs, 
and its sensitivity using simulated nasopharyngeal specimens.’ We 
argue that the ELISA may be used in animal field studies. In the present 
study, results of the N protein ELISA showed a concordance with those 
obtained from nucleic acid detection using dromedary nasal samples, 
indicating that the ELISA was highly specific. It is notable that there is 
one amino acid difference (L178V) between the N protein sequence of 
the MERS-CoV isolate of the present study and the original MERS- 
CoV isolate we used to generate the monoclonal antibodies for the N 
protein ELISA. This N protein ELISA would probably be useful in 
laboratories in the Middle East and North Africa, where resources 
and expertise may be limited. Further refinements of the test and 
the development of a rapid antigen detection kit for rapid diagnoses 
are warranted. 
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